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Abstract 

Centrality, rapidity, and transverse momentum dependence of hadron elliptic flow is 
studied in Au+Au collisions at BNL RHIC energies within the microscopic quark- 
gluon string model. The QGSM predictions coincide well with the experimental data 
at y/s = 130 A GeV. Further investigations reveal that multi-Pomeron exchanges and 
hard gluon-gluon scattering in primary collisions, accompanied by the rescattering 
of hadrons in spatially anisotropic system, are the key processes needed for an 
adequate description of the data. These processes become essentially important for 
heavy-ion collisions at full RHIC energy y/s = 200A GeV. 
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One of the main goals of experiments on heavy-ion accelerators, including 
the Relativistic Heavy-Ion Collider (RHIC) at BNL, which is operating since 
1999, and the coming Large Hadron Collider (LHC) at CERN, is to deter- 
mine the equation of state (EOS) of nuclear matter under extreme conditions 
and to probe the deconfmement phase transition to the new state of matter 
dubbed quark-gluon plasma (QGP) (see, e.g., [1]). It has been understood 
long ago that nuclear matter expanding in the direction perpendicular to the 
beam axis of the colliding nuclei can carry information about the early stage 
of the reaction and subsequent evolution of the system [2]. This phenomenon 
known as "collective flow" has been intensively studied both theoretically and 
experimentally during the last two decades (see [3-5] and references therein). 
Nowadays the expansion in Fourier series is usually applied to study the az- 
imuthal distribution of particles [6,7]: 
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where <p is the azimuthal angle, and pt and y is the transverse momentum 
and the rapidity of a particle, respectively. The first term in square brackets 
represents the isotropic radial flow, while the others are referred to anisotropic 
flow. The first Fourier coefficient Vi = (cos0) = (p x /pt) is called directed 
flow, and the second one v 2 = cos (20)) = {(p x /pt) 2 — (Py/Pt) 2 ) is dubbed 
elliptic flow. Here the beam axis is labelled as z-axis. Together with the impact 
parameter x-axis, it defines the reaction plane perpendicular to the y-axis. To 
study centrality, (pseudo)rapidity, and transverse momentum dependence of 
the anisotropic flow the following triple differential distributions integrated 
over any two of variables in given intervals are useful 
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'f , , / ) d 3 N , 

v n {Xi) = J cos(n0J- — dxjjti J J — — dxjjti , (2) 



r (i) 



where the variables Xj, j = 1,2,3 denote the impact parameter 6, transverse 
momentum p t , and (pseudo)rapidity {rj)y. 

The importance of elliptic flow as a probe of the hot and dense nuclear phase 
at high energies has been discussed first in Ref. [8]. Both microscopic and 
macroscopic calculations seem to support the idea that elliptic flow is devel- 
oped at the very early stage of a nuclear collision [9-15]. Several features, 
such as a characteristic "kinky" structure in the excitation function of elliptic 
flow [16], have been proposed to search for the QGP phase. The first data on 
elliptic flow of charged particles in the midrapidity range of Au+Au collisions 
at v^s = 130^4 GeV, measured by the STAR Collaboration, became available 
recently [17]. It appears that the microscopic transport cascade models, e.g. 
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ultra-relativistic quantum molecular dynamics (UrQMD) [18] or RQMD [16], 
predict too weak values of V2(b, rj) [19,20], while hydrodynamic models over- 
predict the measured elliptic flow by about 20-50 % [17]. In hydrodynamics 
the EOS can be softened by the introduction of a QGP phase and hadron 
resonance-rich matter [21,22]. In the present paper we demonstrate that the 
elliptic flow at RHIC can quantitatively be described within cascade string 
models by the mechanism of string excitations due to colour exchange. This 
procedure is different from the FRITIOF model of longitudinal excitations 
of strings [23] which is the basic part of the string sector in several ultra- 
relativistic transport models. For the further analysis the quark-gluon string 
model (QGSM) [24] is employed. 

In the hadron-hadron (hh) collision part of the QGSM statistical weights, 
hadron structure functions, and leading quark fragmentation functions have 
been obtained in Ref. [25] within the Gribov-Regge theory (GRT) [26]. This 
enables one to choose subprocesses of string excitations, calculate mass and 
momentum of a string, and simulate the string fragmentation into hadrons 
properly. Since a hadron is represented by a vector in the Fock space of con- 
stituents, namely valence and sea quarks, diquarks, gluons, and their antis- 
tates, the strings produced in hadronic collisions originate from different in- 
teractions between the constituents. Scattering of strings is approximated by 
scattering of the valence quarks and diquarks at the end of the strings, and of 
hadrons that are produced in string decays. Similar to the RQMD model the 
valence quarks and diquarks are allowed to interact promptly, while the newly 
produced hadrons can scatter only after a formation time related to the time 
needed to break a string. Together with subprocesses with quark annihilation 
and exchange associated with the Reggeon exchanges in GRT, the model in- 
cludes subprocesses with colour exchange connected to one or more Pomeron 
exchanges in elastic amplitudes in the GRT. The hard gluon-gluon scatter- 
ing with large momentum transfer and the so-called semihard processes with 
quark and gluon scattering are included in the QGSM as well [27]. The inelas- 
tic hh cross section <7j n (s) can be calculated via the real part of the eikonal 
u(s, b) 



(s) = 2tt J {l - exp [-2u R (s, b)] } bdb . (3) 



Here s is the center-of-mass energy of the reaction. The eikonal can be pre- 
sented as a sum of three terms corresponding to soft and hard Pomeron 
exchange, and triple Pomeron exchange, which is responsible for the single 
diffraction process, 

u R (s, b) = uf oft (s, b) + uL(s, b) + uf nple {s, b) . (4) 
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Using the Abramovskii-Gribov-Kancheli (AGK) cutting rules [28] the inelastic 
cross section of hh interaction can be presented as 

(5) 
(6) 

k 

The last equation enables one to determine the number of cut soft and hard 
Pomerons, i.e., the number of strings and hard jets. The single Pomeron ex- 
change, that can be represented by a cylinder-type diagram [25], leads to the 
formation of two quark-diquark or quark-antiquark strings. With rising energy 
the processes with multi-Pomeron exchanges become more and more impor- 
tant. The contribution of the cylinder diagrams to the scattering amplitude 
increases as s a H°) _1 ? while that of the so-called chain diagrams corresponding 
to n-Pomeron exchanges (n > 2) rises as s n ' ap ' ^ 1 ' with a P (0) > 1 being the 
intercept of a Pomeranchuk pole. Diagrams corresponding to hard gluon-gluon 
scattering and double Pomeron exchange are shown as an example in Fig. 1(a) 
and (b), respectively. An overview and recent development of the GRT can be 
found in Ref. [29] . The cascade (rescattering of secondaries) is introduced in 
the QGSM as well. There are several other models which employ the colour 
exchange mechanism for string excitations, e.g., the dual parton model (DPM) 
[30] and the very energetic nuclear scattering model (VENUS) [31]. 

The centrality, transverse momentum, and rapidity dependence of anisotropic 
flow in the QGSM at energies from AGS to SPS has been studied in Refs. [32- 
34]. Here we will focus on study of elliptic flow in minimum bias Au+Au 
collisions at two energies, y/s = 130A GeV and 200A GeV, available at 
RHIC. Rapidity and pseudorapidity distributions of elliptic flow of pions and 
charges particles are depicted in Fig. 2. For both energies elliptic flow displays 
strong in-plane alignment in accordance with the predictions of Ref. [8]. In the 
mid(pseudo)rapidity the flow is almost constant. At y/s = 130A GeV it rises 
up slightly at \y\, \r]\ ~ 1.5, and then drops with the increasing rapidity. The 
mean value of the < 1-3) equals 4.16 ± 0.5 %, that is very close to the 

value 4.5 ± 0.4 % measured by the STAR Collaboration [17]. Pseudorapidity 
dependence of the elliptic flow of charged particles in the whole rj range is also 
in a good agreement with the preliminary results reported by the PHOBOS 
Collaboration [35]. The QGSM predicts that at full RHIC energy elliptic flow 
of charged particles will increase further to < 1-3) — 4.9 ± 0.5%. To 

elaborate on the influence of hard processes and multi-Pomeron exchanges on 
the elliptic flow formation the flow caused by the subprocesses without the 



i,j,k=0;i+j+k>l 
oo 

o'ijk(s) = 2tt J bdbexp —2u R (s,b) 
o 



2 < oft (s,b) 2u R ard (s,b) 2u? riple (s,b) 



4 



hard and multichain contributions is also plotted in Fig. 2. It seems that in 
Au+Au collisions at y/s = 130^4 GeV the magnitude of the signal (except of 
the midrapidity range) can be reproduced without the many-string processes. 
At t/s = 200A GeV their role becomes more significant, because the elliptic 
flow caused by other subprocesses cannot exceed the limit of 3—3.5% . Here 
it is important to stress that the multichain diagrams alone, without rescat- 
tering, cannot affect the elliptic flow at all. The flow increases solely due to 
secondary interactions of produced particles in spatially asymmetric system. 

Figure 3 presents the centrality dependence of elliptic flow of charged particles. 
Since the centrality of events in the experiment has been determined via the 
ratio of charged particle multiplicity to its maximum value N c f l /(N c i l ) max , 
we compare the v 2 [N c h/{N c h) max ] signal with the original impact parameter 
dependence v 2 (b). One can see that the ratio N ch /(N ch ) max is a good criterion 
of the event centrality except of the very central region with b < 2.5 fm, 
where the multiplicity depends weakly on the impact parameter. However, as 
a function of the impact parameter b elliptic flow is saturated at b w 8 fm for 
both energies, while as a function of the multiplicity ratio it increases nearly 
linearly with decreasing multiplicity up to N ch /(N ch ) max ps 0.2. As expected, 
the flow in the midrapidity region is caused mainly by pions. The magnitude 
of the pionic flow in the QGSM calculations is twice as large as obtained, e.g. 
with RQMD [20]. Without the many-string processes the QGSM is able to 
describe the flow only in central and semicentral collisions. It predicts a drop 
of the elliptic flow as the reaction becomes more peripheral, which is similar 
to the predictions of other string models [19,20]. 

The difference in centrality dependences of elliptic flow in calculations with 
and without multi-string processes can be explained as follows Q Hard jets 
and many-Pomeron exchanges produce more particles, thus giving rise to more 
rescatterings that drive the system toward thermal and chemical equilibrium. 
The centrality dependence of elliptic flow in this case is very close to that ob- 
tained at the hydrodynamic limit first in Ref. [8]. In contrast, the less number 
of produced secondaries would bring the system closer to the low density limit 
[14], which means incomplete (if any) thermalization of the system. In this 
regime elliptic flow in peripheral collisions vanishes much earlier compared to 
the hydrodynamic one [14,36]. 

The transverse momentum dependence of elliptic flow is shown in Fig. 4. The 
agreement with the experimental data is good, although the model slightly 
overpredicts the elliptic flow of charged particles at p t < 0.5 GeV/c. Also, the 
dip in the excitation function of elliptic flow at pt ~ 1 GeV/c looks peculiar. 
Nevertheless, both effects are attributed to the lack of heavy resonances in the 
model. This, however, cannot affect the rapidity and centrality distributions 
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where the spectra of particles are integrated over the whole transverse mo- 
mentum range. Similar to the centrality dependence, the calculations without 
the many-string processes can explain the experimental v 2 (pt) data only at 
Pt < 0.8 GeV/c. To describe the further increase of elliptic flow with rising p t 
the full set of diagrams is required. At p t > 1.8 GeV/c elliptic flow saturates in 
accord with experimental results [37]. Note that even without the many-string 
processes the elliptic flow of hadrons in the QGSM is stronger than that of 
the string models based on the FRITIOF routine. This is due to the fact that 
the string stretched between the constituents belonging to different hadrons 
is not parallel to the beam axis, thus giving the additional transverse push to 
secondaries. 

In conclusion, the microscopic quark-gluon string model based on the colour 
exchange mechanism is able to reproduce quantitatively experimental data on 
elliptic flow in Au+Au collisions at RHIC energies quite well, i.e. the micro- 
scopic models are still on the market for ultrarelativistic heavy-ion collisions. 
The key processes in the formation of strong elliptic flow at collider ener- 
gies appear to be hard Pomeron exchange, resulting in the creation of hard 
partonic jets, and soft multi-Pomeron exchanges, accompanied by the subse- 
quent rescattering of secondaries in spatially asymmetric system. Note, that 
directed and elliptic flow strongly depends on the mean free path of hadrons, 
which is inversely proportional to the particle density and interaction cross sec- 
tion. The larger the density (or cross section) of particles participating in sec- 
ondary rescattering, the larger the elliptic flow [38]. Therefore, multi-Pomeron 
exchanges and hard processes increase the yield of secondaries at the initial 
stage of the reaction thus enhancing the elliptic flow. Without the many-string 
processes the model cannot describe the rise of elliptic flow V2{b) for peripheral 
events, while the magnitude of the flow vf l {y) in Au+Au collisions at 130,4 
GeV and 200^4 GeV cannot exceed the limit of 3.5%. For the v 2 (pt) distri- 
bution the QGSM predicts the saturation of elliptic flow at p t > 1.8 GeV/c. 
This finding agrees quantitatively with the experimental results and quali- 
tatively with the predictions of other string models, but deviates from the 
hydrodynamic calculations. More comparison is needed to distinguish unam- 
biguously between hydro- and string- approaches, and to find clear signals of 
the quark-gluon plasma formation. 
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Fig. 1. (a) String formation in hard gluon-gluon scattering and soft Pomeron ex- 
change in nucleon-nuclcon (NN) collision, (b) Formation of four strings as a result 
of double Pomeron exchange in NN collision. 
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Fig. 2. Elliptic flow of pions (solid circles) and charged particles (solid squares) as 
a function of rapidity y and pseudorapidity r] in minimum bias Au+Au collisions in 
the QGSM at y/s = 130 A GeV (upper panels) and 200,4 GeV (lower panels). Open 
circles and squares show the elliptic flow in the model version without the hard 
processes and multi-chain diagrams. Hatched area indicates the data measured by 
the STAR Collaboration [17], open triangles denote the preliminary results of the 
PHOBOS Collaboration [35]. 
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Fig. 3. Elliptic flow of charged particles (solid circles) with \rj\ < 1.3 as a function 
of the ratio n c h/(n c h)max (left panels) and the impact parameter b (right panels) 
in minimum bias Au+Au collisions at yfs = 130^4 GeV (upper row) and 200.4 
GeV (lower row). The flow of charged pions with \y\ < 1 in the model is shown 
by open triangles, open circles denote the experimental data from Ref.[17]. Open 
squares indicate the elliptic flow in the model version without the hard processes 
and multi-chain diagrams. 
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Fig. 4. Transverse momentum dependence of the elliptic flow of charged particles 
with | r]\ < 1.3 in minimum bias Au+Au collisions at 130^4 GeV and 200^4 GeV. 
Data from Ref.[17] are shown by open squares. 
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